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Abstract 

We present a perturbative QCD factorization formalism for the production of heavy quarkonia of 
large transverse momentum pT at collider energies, which includes both the leading power (LP) and 
next-to- leading power (NLP) contributions to the cross section in the mq/p'^ expansion for heavy 
quark mass mq. We estimate fragmentation functions in the non-relativistic QCD formalism, and 
reproduce the bulk of the large enhancement found in explicit NLO calculations in the color singlet 
model. Heavy quarkonia produced from NLP channels prefer longitudinal polarization. 
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Introduction — More than thirty-five years after the discovery of the J/ip , the production 
of heavy quarkonia remains a key subject in strong interaction physics [2]. The inclusive 
production of pairs of charm or bottom quarks, with masses niq ^ Aqcd, is an essentially 
perturbative process, while the subsequent evolution of the pair into a physical quarkonium is 
nonperturbative. Different treatments of the transformation from heavy quark pair to bound 
quarkonium are given in various formalisms, most notably, the color singlet model (CSM), 
the color evaporation model (CEM), and non-relativistic QCD (NRQCD . The current status 
of theory and experiment has been summarized very recently in Ref . ^ \ . 

For the NRQCD formalism [s], Q], small, color-octet production matrix elements can 
provide good fits to high px inclusive hadron collider cross sections for J/ip and T, but a 
complete description remains elusive. Polarization in particular remains a challenge, along 
with the surprisingly high rate of associated production in electron-positron annihilation 
[2I. Although plausible arguments for the use of NRQCD at leading power (p^"^) have been 
around since the beginning of the formalism \4 , issues of gauge invariance and infrared 
cancellation are still not completely settled jsl,]^]. Meanwhile, large next-to-leading order 
(NLO) and potentially large next-to-next-to leading order (NNLO) corrections to high-p^ 
cross sections [3] in the CSM have attracted attention. The size of these color singlet 
cross sections seems to upset long- held expectations for gluon fragmentation/color octet 
dominance, and indeed, they appear in diagrams that fall off like p^^ in da/dp"^, compared to 
the leading, p^"^ behavior associated with gluon fragmentation. These developments suggest 
that we must widen the formalism for quarkonium production beyond leading power. 
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FIG. 1. J /ip production through a single parton (a gluon) fragmentation (left) and fragmentation 
of a heavy quark pair (right). 
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In this letter, we emphasize the expansion of the production rate of the heavy quark pairs 
first in the large scale pr and then in the coupling constant as- We present a perturbative 
QCD factorization formalism, accurate to the first nonleading power in ttlq/p^ that incor- 
porates both leading power gluon fragmentation and direct production of heavy quarks at 
short distances with subsequent fragmentation, as illustrated in Fig. [T] jsl- These figures are 
shown in cut diagram notation, in which the amplitude and complex conjugate are com- 
bined into a forward scattering diagram in which the final state is identified by a vertical 
line. When px ^ "^g, the rate to produce a gluon (or in general a single parton) at distance 
scale 1/pt; which fragments into a heavy quark pair and thence into a physical quarkonium, 
is characterized by a behavior. We refer to this as the leading power (LP) contribution. 
The perturbative production of a coUinear heavy quark pair directly at the short distance 
scale 1/pt is suppressed by l/p^ relative to the production rate of a single parton at the 
same pt- We therefore refer to it as a next-to-leading power (NLP) contribution. As we 
will see below, however, the probability for a such heavy quark pair to evolve into a heavy 
quarkonium is naturally enhanced compared to that of a single parton. This can promote 
the NLP channel to phenomenological interest, despite its suppression by two powers oipT- 



The physical heavy quarkonium is likely formed long after the heavy quark pair was pro- 
duced [9| . Both the LP and NLP production channels in Fig. [1] can therefore be considered 
as fragmentation contributions [2]. Similar to factorization at LP the NLP contribution 
to the production rate is factorized into perturbatively calculable short-distance coefficient 
functions for producing the heavy quark pair convoluted with non-perturbative but univer- 
sal long-distance fragmentation functions for the pair to become a heavy quarkonium Q]. 
The universality of the fragmentation functions can in principle be tested when we compare 
the data of heavy quarkonium production from processes with different short-distance co- 
efficients. As a consequence of this perturbative factorization, the short-distance coefficient 
functions at both the LP and NLP capture the dynamics at the distance scale 1/pt, and 
are insensitive to the details of the produced heavy quarkonium states. 



Factorization formalism — For the production of a heavy quarkonium state H of momentum 
p, A[Pi) + B[P2) — 7- H[p) + X, the leading contribution from the channels in Fig. [T]can be 
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summarized in an extended factorization formula given schematically by 

daA+B-,H+xip) ~ (1) 

^daA+B^f+x{Pf =p/z)(g) DH/f{z,mQ) 
f 

+ daA+B^[QQi.)]+x{pi'^±0/2z,p{l±C)/2z) 

[QQ{k)] 

where the factorization scale dependence is suppressed. The first (second) term on the right- 
hand side gives the contribution of LP (NLP) in mq/pT- In the first term, daA+B^f+x{Pf) 
is the cross section to produce an on-shell parton of flavor / at short distances, which 
contains all of the information about the incoming state and includes convolutions with 
parton distributions when A or 5 is a hadron. The sum over / runs over all parton flavors, 
and Dfj/fiz, rriQ) is the fragmentation function for a heavy quarkonium state H from parton 
/ with momentum fraction z B- For H a J/ip or other heavy quarkonium, the dominant 
channel at hadron colliders is gluon fragmentation, f = g. The second term on the right-hand 
side in Eq. ([1]) is suppressed by relative to the first, and the quark-pair fragmentation 
function, (, mq), has units of mass squared, which are compensated by large 

invariants from the hard-scattering function, daA+B~>[QQ{K)]+x{p{^ i 0/2^,^(1 ± (')/2z), 
which describes production of an on-shell, coUinear heavy quark pair. The momentum 
fractions z, ( and (' are defined as 

+ +1+C + +1-C 
" ^ ^ ' Pq=p ' 

^« ^ ^ 

By analogy to the single-parton case, z measures the fractional momentum of the pair carried 
by the observed quarkonium in this leading region, which is the same on both sides of the 
cut in Fig. [TJ Parameters C and C' characterize the sharing of the pair's momentum between 
the heavy quark and antiquark on either side of the cut in the figure. In principle, these need 
not be the same. The ® in Eq. (HI) represents the convolution over the partons' momentum 

n 

fractions [10]. 

In Eq. ([T]), the perturbative hard parts at both the LP and NLP capture the QCD 
dynamics at distance scale of l/p^, are independent of heavy quark mass, and are the same 
for the production of all heavy quarkonium states. For the LP, the hard parts are effectively 



the same as those for inclusive production of any single hadron at high pt and are available 



to both leading order (LO) and next-to- leading order (NLO) in 

For the hard parts at NLP, we only need to calculate the rate to produce a heavy quark 
pair with zero relative transverse momentum, since the effect of relative transverse mo- 
mentum will be further suppressed in l/pr- The hard parts can be perturbatively calcu- 
lated order-by-order in powers of by applying the factorization formula in Eq. ([T]) to 
the production of an asymptotic state of a heavy quark pair [QQ{k)] of momentum p with 
various spin and color quantum numbers. For example, the hard part for the subprocess, 
q + q —> [QQ{k)] + g, can be derived by applying Eq. ([1]) to heavy pair production in 
quark- ant iquark scattering, 

-dal^ig^g^giPg = p/z) ® V^^Q^^^^^p, tuq) , 

where the superscript "(n)" indicates the order in as- The first term on the right of Eq. ([3]) 
is the differential cross section for the subprocess given by the diagrams in Fig. |2|^a) and (b) 
plus three more diagrams. When px ^ mg, the square of diagram (a) clearly contributes 
to the NLP, while the the square of diagram in (b) contributes to both LP and NLP. We 
note that the interference between the two contributes only to NLP. The second term in 
Eq. (|3]) exactly removes the LP piece of this subprocess, where dag^g_^g_^g{pg) is the LO 
differential cross section for g + g — )■ g{pg) +g and V^^Q^^.^y^{z,mQ) is the LO fragmentation 
function given by the (cut) diagram in Fig. Et^c). Since the second term on the right of 
Eq. ([3]) removes a power mass singularity ~ I/itlq, it is important to keep the heavy quark 
mass when evaluating partonic diagrams and to set mq — )■ (only) after carrying out the 
subtraction jlo| . 

The predictive power of the factorization formula in Eq. ([1]) relies on the universality of the 
fragmentation functions. The single parton fragmentation functions Dn/fiz, mq) are defined 
in Ref. jsl. The operator definition for ^^///[qq(k)] {z, (, mq) depends on k, which represents 
the pair's color and spin. For the heavy quark pair moving in the direction with light- 

cone momentum components, p^ = (p"*", {2mqy / 2p'^ , 0±), there are singlet (1) and octet (8) 
color states, and four spin states described by relativistic Dirac spin projection operators: 
7+75/4^+, 7+/4p"'', and 7+7Y4p+ with i = 1,2, for effective axial vector (a), vector (v), and 
tensor (t) "currents" , respectively [lOi] . As an example, the operator definition of the axial 







(b) (c) 
FIG. 2. Sample Feynman diagrams for g + g — )• [QQ{k)]{p) + g subprocess. The diagram in (a) 
contributes to the short-distance coefficient at NLP while the diagram in (b) contributes to both 
the LP and NLP hard parts. The cut diagram in (c) gives the LO gluon to heavy quark pair 
fragmentation function. 



vector/octet fragmentation function can be written as 
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where we have suppressed dependence on a factorization scale. We also suppress gauge links 
along the minus light cone, inserted between repeated color indices, which provide a gauge 
invariant definition of the operator [l^. The links are in adjoint representation for index 
a, and fundamental representation for i,j,l,k. (We note that the color matrices may 
be taken at any points along the light cone.) Overall, the pair fragmentation functions are 
given by matrix elements of nonlocal operators, and the form is very similar to operator 
definitions of single parton fragmentation functions, simply replacing the parton field by the 
product of quark fields. They are also reminiscent of hadronic wave functions that connect 
multiple partons to the hard scattering in the factorized expressions for elastic amplitudes 

Q. 

Similar to the single parton fragmentation functions, heavy quark pair fragmentation 
functions like the one defined in Eq. (jl]) are nonperturbative, but, universal. They are boost 
invariant, and require renormalization. They thus evolve in the usual sense, and depend 



upon a factorization scale, chosen to match the short-distance scale of the problem. It is 
natural to think of the choice of factorization scale as the same for NLP as for LP, and in 
general, a change of scale could mix them. This leads to a closed set of general evolution 
equations for both LP and NLP fragmentation functions, which will be addressed elsewhere 
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We note that there are further additive corrections at the power of such as those in- 
volving twist-4 parton distributions, or twist-4 light-parton fragmentation functions. These 
corrections, however, can be considered small since they do not introduce natural enhance- 
ments that scale with rriQ in the probability to form a quarkonium, to compensate the 
suppression of l/p"^. 

Cross section and polarization — It is the fragmentation functions that determine the ab- 
solute normalization of perturbative calculations using the factorization formula in Eq. ([1]), 
as well as the differences in the production rate between various quarkonium states and 
their polarizations. Since only pairs with small relative momentum are likely to form bound 
quarkonia, we may apply the basic NRQCD factorization hypothesis for heavy quarkonium 
production to these fragmentation functions to reduce the unknown functions to a few uni- 
versal constants in the form of NRQCD matrix elements, 

c 

c 

where the functions d are calculable, is a factorization scale, and (0|^q(c)]) local 



NRQCD matrix elements {4]. Although we cannot provide a full proof for the NRQCD 
factorization in Eq. ([5]), it is reasonable to evaluate the coefficient functions in Eq. to 
estimate the properties of fragmentation functions 5|. 

We can use this formalism to help understand the source of the surprisingly large correc- 
tions to J/tp production in the CSM at NLO and NNLO with predominantly longitudinal 



polarization, in contrast to the small and transversely polarized LO 13|, [ij]. To do so, 
we use Eq. ([T]) to calculate NLP J/ip cross sections at LO in from a color singlet '^Si 
pair. We estimate the 0{as) fragmentation functions from the diagrams in Fig. [3l where 
the upper lines are fixed at p/2 and Dirac indices are contracted with an NRQCD singlet 
projection with matrix element essentially equivalent to the CSM. At order of a^, only the 



[QQ{a8)] state fragments into a color singlet ^Si heavy quark pair, while both \QQ(y8)] 
and [QQ(t8)] states give vanishing contributions due to charge conjugation symmetry |lO|. 
Defining r{z) = z^/i^/(4m^(l — z)"^), we have 



1 (Offs.)) 
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(6) 



where {T^'^) is the pair-fragmentation function with the pair longitudinally (transversely) 
polarized, and 



(7) 



with a, 6 



-1,1, and total unpolarized contribution, V 



2V 



T _ 

[QQ(a8)]^J/V 



+ 



■^^Q(a8)]->j/t/'' deriving Eq. ([H]), we renormalized the UV divergence by a cutoff /i on 
transverse momenta. Quark masses make an IR cutoff on transverse momenta unnecessary. 
Other renormalization schemes give similar results and will be discussed elsewhere 10 1. 



With these fragmentation functions, and LO hard scattering functions Q0|, we can compute 



NLP cross sections that can be compared to NLO GSM cross sections [7|, \mm. Using 
the CTEQ6L parton distributions [17], fi = px for factorization and renormalization scales, 
rric = 1.5 GeV, and (O^sg^^) = 1.32 GeV^, we calculate and plot the unpolarized cross 
section as a function pj- in Fig. H] (upper panel). We see that the factorization-based cross 
section shows much the same enhancement above the LO GSM as the full NLO result in 
GSM. Actually, our calculated NLP cross section with the estimated 0{as) fragmentation 
functions reproduces as much as eighty percent of the full NLO GSM result. The difference 
should be due mainly to 0{as) corrections to the hard scattering functions from the color 
singlet [QQ(^l)] channel. With higher order fragmentation functions from other [QQ{k)] 
states, the full NLP cross section at LO in could be larger than the solid line in Fig. SI 
Similarly, we calculate the J/ip polarization (lower panel), as measured by the parameter 



a 



a — a )/{a + a ) in terms of transverse (longitudinal) cross section a {a ). Our 
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FIG. 3. Leading order Feynman diagrams represent the fragmentation of a heavy quark pair to 
another heavy quark pair. 
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FIG. 4. J/-!/; cross section (upper panel) and polarization (lower panel) as a function of px- The 
solid lines are computed from the NLP term in Eq. ([1]) with fragmentation functions given by Eq. 
([6]); the dashed line is LO CSM cross section as in Ref. 13], and all in [nb/GeV]. 

result is consistent with that in Refs. We regard these results as compelling evidence 

for the phenomenological relevance of the power expansion. 

Summary — We have argued that a practical strategy for the phenomenology of heavy 
quarkonium production at high pt is to expand the cross sections first in the large scale 



Pt, and then in the couphng, ag, and have presented a new perturbative QCD factorization 
formahsm for quarkonium production including the first non-leading powers in mg/pT- This 
approach enables us to resum perturbative logarithms into the fragmentation functions, to 
analyze systematically the influence of the larger color singlet matrix elements despite their 
suppressed pT-dependence, and to resolve some of the mystery associated with the discovery 
of large high order corrections to color singlet cross sections. We have found that heavy 
quarkonia produced from pair fragmentations are likely to be longitudinally polarized, in 
contrast to single parton fragmentation. The observed quarkonium polarization should be 
a consequence of the competition of these two leading production channels. 
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